Equation of state of atomic systems beyond s-wave determined by the lowest order 
constrained variational method: Large scattering length limit 
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Dilute Fermi systems with large s-wave scattering length a s exhibit universal properties if the 
interparticle spacing r greatly exceeds the range of the underlying two-body interaction potential. 
In this regime, r is the only relevant length scale and observables such as the energy per particle 
depend only on r (or, equivalently, the energy Efg of the free Fermi gas). This paper investigates 
Bose and Fermi systems with non- vanishing angular momentum I using the lowest order constrained 
variational method. We focus on the regime where the generalized scattering length becomes large 
and determine the relevant length scales. For Bose gases with large generalized scattering lengths, 
we obtain simple expressions for the energy per particle in terms of a /-dependent length scale 
which depends on the range of the underlying two-body potential and the average interparticle 
spacing. We discuss possible implications for dilute two-component Fermi systems with finite I. 
Furthermore, we determine the equation of state of liquid and gaseous bosonic helium. 

PACS numbers: 
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I. INTRODUCTION 

The experimental realization of dilute degenerate Bose 
and Fermi gases has led to an explosion of activities in the 
field of cold atom gases. A particularly intriguing feature 
of atomic Bose and Fermi gases is that their interaction 
strengths can be tuned experimentally through the ap- 
plication of an external magnetic field in the vicinity of 
a Feshbach resonance [l|, 0]. This external knob allows 
dilute systems with essentially any interaction strength, 
including infinitely strongly attractive and repulsive in- 
teractions, to be realized. Feshbach resonances have been 
experimentally observed for s-, p- and d-wave interacting 
gases [1, H], [H, @, 0] and have been predicted to exist also 
for higher partial waves. 

A Feshbach resonance arises due to the coupling of two 
Born-Oppenheimer potential curves coupled through a 
hyperfinc Hamiltonian, and requires, in general, a multi- 
channel description. For s-wave interacting systems, Fes- 
hbach resonances can be classified as broad or narrow [1] . 
Whether a resonance is broad or narrow depends on 
whether the energy width of the resonance is large or 
small compared to the characteristic energy scale, such 
as the Fermi energy or the harmonic oscillator energy, 
of the system. In contrast to s-wave resonances, higher 
partial wave resonances are necessarily narrow due to the 
presence of the angular momentum barrier Q. This pa- 
per uses an effective single channel description to investi- 
gate the behaviors of strongly-interacting Bose and Fermi 
systems with different orbital angular momenta. 

In dilute homogeneous Bose and Fermi gases with large 
s-wave scattering length a s , a regime has been identi- 
fied in which the energy per particle takes on a universal 
value which is set by a single length scale, the average 
ToL 111 . H2I . In this so-called uni- 



interparticle spacing r a 
tary regime, the length scales of the s-wave interacting 
system separate according to \a s \ ^> r a ^> R, where R 
denotes the range of the two-body potential. The en- 



ergy per particle Eb,o/N (the subscripts ll B" and "0" 
stand respectively for "boson" and "s-wave interacting" ) 
for a homogeneous one-component gas of bosons with 
mass m in the unitary regime has been calculated to be 
E B ,a/N sa 13.3 ft 2 nf/ to using the lowest order con- 
strained variational (LOCV) method [l^|. The energy 
Eb,o/N at unitarity is thus independent of a s and R, 
and depends on the single length scale r a through the 
boson number density Ub, t = (47rrts/3) -1 / 3 . How- 
ever, Bose gases in the large scattering length limit are 
expected to be unstable due to three-body recombina- 
tion [H, 0, [H [g. 

On the other hand, the Fermi pressure prevents the 
collapse of two-component Fermi gases with equal masses 
and equal number of "spin-up" and "spin-down" fermions 
with large interspecies s-wave scattering length [l(| EL 
[l7l [H[ . At unitarity, the energy per particle is given by 
E F ,o/N w QA2E FG , where E FG = (3 / 10) (h 2 k F /m) de- 
notes the e nergy p er particle of the non-interacting Fermi 
gas (T9I I20I I21I |22| . The Fermi wave vector k F is related 
to the number density of the Fermi gas by n F = kp/3ir 2 , 
which implies that E Fj q/N depends on r a but is indepen- 
dent of a s and R. We note that the inequality \a s \ 3> r a 
is equivalent to l/(k F \a s \) <C 1. 

This paper investigates Bose and Fermi systems with 
large generalized scattering lengths using the LOCV 
method. For p- and d-wave interacting Bose systems, 
we define the unitary regime [23| through the inequali- 
ties \ai(E re i) \ ^> & 3> R, where £j denotes a ^-dependent 
length scale given by the geometric combination of r and 
R, i.e., = ro 1 '^i?'/ 4 , and E re i the relative scatter- 
ing energy. The generalized energy-dependent scattering 
length ai(E re i) 0, [H [2^] characterizes the scattering 
strength (see below) . We find that the energy of p-wave 
interacting two-component Bose gases and d-wave inter- 
acting one- and two-component Bose gases at unitary is 
determined by the combined length £1 . While Bose gases 
with higher angular momentum in the unitary regime are 
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of theoretical interest, they are, like their s-wave cousin, 
expected to be unstable. We comment that the energet- 
ics of two-component Fermi gases with large generalized 
scattering length may depend on the same length scales. 

Furthermore, we consider s-wave interacting Bose sys- 
tems over a wide range of densities. Motivated by two re- 
cent studies by Gao [27], HH , we determine the energy per 
particle Eg /N of the Bose system characterized by two 
atomic physics parameters, the s-wave scattering lengh 
a s and the van der Waals coefficient C§. Our results lead 
to a phase diagram of liquid helium in the low-density 
regime that differs from that proposed in Ref. [28| . 

Section |TT] describes the systems under study and in- 
troduces the LOCV method. Section QTI] describes our 
results for dilute s-wave interacting Bose and Fermi sys- 
tems and for liquid helium. Section IIVI considers Bose 
and Fermi systems interacting through i-wave (I > 0) 
scattering. Finally, Section fVl concludes. 



II. LOCV METHOD FOR BOSONS AND 
FERMIONS 

This section introduces the three-dimensional Bose 
and Fermi s yst ems under study and reviews the LOCV 
method Hfi EH, H2. The idea of the LOCV method 
is to explicitly treat two-body correlations, but to ne- 
glect three- and higher-body correlations. This allows 
the many-body problem to be reduced to solving an 
effective two-body equation with properly chosen con- 
straints. Imposing these constraints makes the method 
non- variational, i.e., the resulting energy does not place 
an upper bound on the many-body energy. The LOCV 
method is expected to capture some of the key physics of 
dilute Bose and Fermi systems. 

The Hamiltonian Hg for a homogeneous system con- 
sisting of identical mass m bosons is given by 



H R = — 



2m 



(i) 



i<j 



where the spherically symmetric interaction potential v 



depends on the relative distance fj 



|rj-r,- . Here, 



Vi denotes the position vector of the ith boson. The 
Hamiltonian Hp for a two-component Fermi system with 
equal masses and identical spin population is given by 



H, 



h 2 h 2 

= -— y v 2 -— y v 2 



(2) 



where the unprimed subscripts label spin-up and the 
primed subscripts spin-down fermions. Throughout, we 
take like fermions to be non-interacting. Our primary 
interest in this paper is in the description of systems 
for which many-body observables are insensitive to the 
short-range behavior of the atom-atom potential v(r). 
This motivates us to consider two simple model poten- 
tials: an attractive square well potential v sw with depth 



V (V > 0), 



,(r) 



-V 




for r < R 
for r > B : 



(3) 



and an attractive van der Waals potential v vc i w with hard- 
core rv, 



v v dw(r) 



-C 6 /r 



for 
for 



r <r c 
r > r c 



(4) 



In all applications, we choose the hardcore r c so that the 
inequality r c <C /3 6 , where (3% = (mCe/H 2 ) 1 ^ 4 , is satisfied. 
The natural length scale of the square well potential is 
given by the range R and that of the van der Waals po- 
tential by the van der Waals length (3q. The solutions to 
the two-body Schrodinger equation for v sw are given in 
terms of spherical Bessel and Neumann functions (impos- 
ing the proper continuity conditions of the wave function 
and its derivations at those r values where the poten- 
tial exhibits a discontinuity) , and those for v v d w in terms 
of convergent infinite series of spherical Bessel and Neu- 
mann functions [33j |. 

The interaction strength of the short-range square well 
potential can be characterized by the generalized energy- 
dependent scattering lengths ai(k), 



ai(k) = sgn[-tan^(fc)] 



tan5;(fc) 



k 2l+1 



l/(2/+l) 



(5) 



where Si(k) denotes the phase shift of the Ith partial 
wave calculated at the relative scattering energy Erei ? 
k = mE re i / h 2 . This definition ensures that ai(k) ap- 
proaches a constant as k — > [34l l35j . For the van der 
Waals potential v v d w , the threshold behavior changes for 
higher partial waves and the definition of a/(fc) has to be 
modified accordingly [34], l35j . In general, for a poten- 
tial that falls off as — r~ n at large interparticle distances, 
ai(k) is defined by Eq. ([5|) if 21 < n — 3 and by 



ai(k) = sgn[- tan 



tan Si(k) 



l/(n-2) 



(6) 



if 2Z > n — 3. For our van der Waals potential, n is equal 
to 6 and ai(k) is given by Eq. for I < I and by Eq. © 
for I > 2. The zero-energy generalized scattering lengths 
ai can now be defined readily through 



ai = lim ar(k). 

fc^O 



(7) 



We note that a new two-body Z-wave bound state ap- 
pears at threshold when |a/| — > oo. The unitary regime 
for higher partial waves discussed in Sec. IIVI is thus, as 
in the s-wave case, closely related to the physics of ex- 
tremely weakly-bound atom-pairs. To uncover the key 
behaviors at unitarity, we assume in the following that 
the many-body system under study is interacting through 
a single partial wave I. While this may not be exactly 
realized in an experiment, this situation may be approx- 
imated by utilizing Feshbach resonances. 
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We now outline how the energy per particle Eb,i/N 
of a one-component Bose system with /-wave interac- 
tio ns |36l| is calculated by the LOCV method [1^, [3(| 
l3lL |32|| . The boson wave function is taken to be a 
product of pair functions fi, 



,r N ) 



(8) 



and the energy expectation value of Hb, Eq. {TJ, is calcu- 
lated using ^b- If terms depending on the coordinates of 
three or more different particles are neglected, the result- 
ing energy is given by the two-body term in the cluster 
expansion, 



E 



B.I 



N 



2 



fl(r) 



V 2 +v(r) Mr) d 3 r. (9) 



The idea of the LOCV method is now to introduce a heal- 
ing distance d beyond which the pair correlation function 
fi is constant, 



Mr >d) = l. 



(10) 



To ensure that the derivative of /; is continuous at r = d, 
an additional constraint is introduced, 



//(r = d) = Q. 



(11) 



Introducing a constant average field A; and varying with 
respect to /; while using that /; is constant for r > d, 
gives the Schrodinger-like two-body equation for r < d, 



-V 2 +v(r) (,•/,(,•)) = .W,(r). 
m 



Finally, the condition 



n B 



fi(r) d 3 r = l 



(12) 



(13) 



enforces that the average number of particles within d 
equals 1. Using Eqs. (flQ]) and (fl2|). the energy per 
particle becomes, 



E 



B.I 



N 



A, 
2 



n B 



v(r)d 3 r. 



(14) 



The second term on the right hand side of Eq. (fl4|) is 
identically zero for the square well potential v sw but con- 
tributes a so-called tail or mean-field energy for the van 
der Waals potential v v d w [13, HI]. We determine the 
three unknown ns, Xi and d by simultaneously solving 
Eqs. (TT2"|) and (fT5|) subject to the boundary condition 
given by Eq. (jlip . Note that ns and d depend, just as 
fi and A;, on the angular momentum; the subscript has 
been dropped, however, for notational convenience. 

In addition to one-component Bose systems, Sec. IIVI 
considers two-component Bose systems, characterized by 
/-wave interspecies and vanishing intraspecies interac- 
tions. The Hamiltonian for the two-component Bose sys- 
tem is given by Eq. with the sum of the two-body 



interactions restricted to unlike bosons. Correspondingly, 
the product wave function is written as a product of 
pair functions, including only correlations between unlike 
bosons. The LOCV equations are then given by Eqs. (TIU|) 
through (Ti"3")) with hb in Eq. (fi"3"|) replaced by Ub/2. 

Next, we discuss how to determine the energy Ep,i/N 
per particle for a two-component Fermi system within 
the LOCV method [H, [H Ml, II . The wavefunction is 
taken to be 



* F (n, . . . ,rv, . . . ) = $ FG Y[ fi(nr), 



(15) 



where $ fg denotes the ground state wavefunction of the 
non-interacting Fermi gas. The product of pair functions 
fi accounts for the correlations between unlike fermions. 
In accord with our assumption that like fermions are non- 
interacting, Eq. (|15[) treats like fermion pairs as uncor- 
rected. Neglecting exchange effects, the derivation of 
the LOCV equations parallels that outlined above for the 
bosons. The boundary conditions, given by Eqs. (|10[) and 
(fTTj) . and the Schrodinger-like differential equation for A;, 
Eq. (fT2|) . are unchanged. The "normalization condition," 
however, becomes 



rip 
~2~ 



f /f(r)d 3 r = 1, 
Jo 



(16) 



where the left-hand side is the number of fermion pairs 
within d. The fermion energy per particle is then the sum 
of the one-particle contribution from the non-interacting 
Fermi gas and the pair correlation energy A; [2C| , 



Ef,i/N — Efg + 



(17) 



This equation excludes the contribution from the tail of 
the potential, i.e., the term analogous to the second term 
on the right hand side of Eq. (fT4|) , since this term is neg- 
ligible for the fermion densities considered in this paper. 

The LOCV solutions for fi, A/ and d for the ho- 
mogeneous one-component Bose system and the two- 
component Fermi system are formally identical if the bo- 
son density is chosen to equal half the fermion density, 
i.e., if ns — tif/2. This relation can be understood by 
realizing that any given fermion (e.g., a spin- up particle) 
interacts with only half of the total number of fermions 
(e.g., all the spin-down fermions). Consequently, the two- 
component Fermi system appears twice as dense as the 
one-component Bose system. The fact that the LOCV 
solutions for bosons can be converted to LOCV solutions 
for fermions suggests that some physics of the bosonic 
system can be understood in terms of the fermionic sys- 
tem and vice versa. In fact, it has been shown previ- 
ously [2fJ that the LOCV energy for the first excited 
gas-like state of s-wave interacting fermions at unitarity 
can be derived from the LOCV energy of the energet- 
ically lowest-lying gas-like branch of s-wave interacting 
bosons [l3|. Here, we extend this analysis and show that 
the ground state energy of the Fermi gas at unitarity can 
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FIG. 1: Energy per particle Eb,o/N as a function of the 
density ns, both plotted as dimensionless quantities, for a 
one-component Bose system interacting through the van der 
Waals potential with s-wave scattering length a 3 = 16.9/3e. 
The dotted line shows the gas branch and the dashed line the 
liquid branch. The minimum of the liquid branch is discussed 
in reference to liquid 4 He in the text. 



be derived from the energetically highest-lying liquid-like 
branch of the Bose system. Furthermore, we extend this 
analysis to higher angular momentum scattering. 



III. s-WAVE INTERACTING BOSE AND 
FERMI SYSTEMS 

Figure [TJ shows the energy per particle Eb,o/N, 
Eq. (TH|) , obtained by solving the LOCV equations for 
a one-component Bose system interacting through the 
van der Waals potential with s-wave scattering length 
a s = 16.9/36- The dotted line in Fig.[TJhas positive energy 
and increases with increasing density; it describes the 
energetically lowest-lying "gas branch" for the Bose sys- 
tem with a s — 16.9/?6 and corresponds to the metastable 
gaseous condensate studied experimentally. The dashed 
line in Fig. [1] has negative energy at small densities, de- 
creases with increasing density, and then exhibits a mini- 
mum; this dashed line describes the energetically highest- 
lying "liquid branch" for a Bose system with a s — 16.9(3q. 
Within the LOCV framework, these two branches arise 
because the Schrodinger-like equation, Eq. (Tl2|) . permits 
for a given interaction potential solutions /o with differ- 
ing number of nodes, which in turn give rise to a host of 
liquid and gas branches [27], [28| . Throughout this work 
we only consider the energetically highest-lying liquid 
branch with n nodes and the energetically lowest-lying 
gas branch with n + 1 nodes. To obtain Fig. [TJ we con- 
sider a class of two-body potentials with fixed a s //?6, and 
decrease the value of the ratio r c //3 6 till Eb,o/N, Eq. (|T4"J) . 
no longer changes over the density range of interest, i.e., 
the number of nodes n of the energetically highest-lying 
liquid branch is increased till convergence is reached. 



FIG. 2: Energy per particle Eb,o/N as a function of the den- 
sity tib for a one-component Bose system interacting through 
the van der Waals potential with s-wave scattering lengths 
a s = 16.9/36 (open circles) and a s = 169/36 (filled circles). 
To guide the eye, dashed and dotted lines connect the data 
points of the liquid and gas branches, respectively. The liq- 
uid branches go to Edimer as the density goes to zero. The 
solid lines show Eb,o/N at unitarity; see text for discussion. 
Compared to Fig. [TJ the energy and density scales are greatly 
enlarged. 

In Fig. [TJ the two-body van der Waals potential is cho- 
sen so that the scattering length of a s = 16.9/36 coin- 
cides with that of the 4 He pair potential [37[ • The liquid 
branch in Fig. [TJ can hence be applied to liquid 4 He, and 
has previously been considered in Refs. [27, 28]. The min- 
imum of the liquid branch at a density of ub = 2.83/3^~ 3 , 
or 1.82 x 10 22 cm -3 , agrees quite well with the experi- 
mental value of 2.18 x 10 22 cm -3 (38|. The corresponding 
energy per particle of —6.56 K deviates by 8.5 % from 
the experimental value of —7.17 K [38j |. This shows that 
the LOCV framework provides a fair description of the 
strongly interacting liquid 4 He system, which is charac- 
terized by interparticle spacings comparable to the range 
of the potential. This is somewhat remarkable consid- 
ering that the LOCV method includes only pair corre- 
lations and that the van der Waals potential used here 
contains only two parameters. 

Open circles connected by a dashed line in Fig. [2] show 
the liquid branch for a s = 16.9/3e in the small density 
region. As the density goes to zero, the energy per par- 
ticle Eb,q/N does not terminate at zero but, instead, 
goes to Edimer/'Z, where Edimer denotes the energy of 
the most weakly-bound s-wave molecule of v v d w ■ In this 
small density limit, the liquid branch describes a gas of 
weakly-bound molecules, in which the interparticle spac- 
ing between the molecules greatly exceeds the size of the 
molecules, and Edimer is to a very good approximation 
given by — fi. 2 /(ma 2 ). As seen in Fig. [51 we find solu- 
tions in the whole density range considered. In contrast 
to our findings, Ref. [28( reports that the LOCV solu- 
tions of the liquid branch disappear at densities smaller 
than a scattering length dependent critical density, i.e., 
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FIG. 3: Scaled interparticle spacing r / (3§ as a function of the 
scaled density ubPs f° r the g as branch of a one-component 
Bose system interacting through the van der Waals poten- 
tial with a a — 169/3e- The horizontal lines show the scaled 
s-wave scattering length a 3 — 169/36 and the range of the van 
der Waals potential, which is one in scaled units (almost in- 
distinguishable from the x-axis). This graph shows that the 
unitary inequalities a s 3> r 3> P& hold for jib larger than 
about lO" 5 /^ 3 . 
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FIG. 4: Scaled energy per particle (Ef,o/N)/Efg as a func- 
tion of \/(kFa s ) for a two-component s-wave Fermi gas inter- 
acting through the square well potential for uf = 
The combined dashed and dash-dotted curve corresponds to 
the BEC-BCS crossover curve and the dotted curve corre- 
sponds to the first excited state of the Fermi gas. The dashed 
and dotted linestyles are chosen to emphasize the connection 
to the gas and liquid branches of the Bose system in Figs. [2] 
and (see text for more details). 



at a critical density of 8.68 x 10~ 7 /3g~ 3 for a s — 16.9/3e- 
Thus we are not able to reproduce the liquid-gas phase 
diagram proposed in Fig. 2 of Ref. [HI, which depends 
on this termination of the liquid branch. We note that 
the liquid branch is, as indicated by its imaginary speed 
of sound, dynamically unstable at sufficiently small den- 
sities. The liquid of weakly-bound bosonic molecules 
discussed here can, as we show below, be related to 
weakly-bound molecules on the BEC side of the BEC- 
BCS crossover curve for two-component Fermi gases. 

We now discuss the gas branch in more detail. Open 
and filled circles connected by dotted lines in Fig. [2] show 
the energy per particle for a s — 16.9/36 and 169 Pe, respec- 
tively. These curves can be applied, e.g., to 85 Rb, whose 
scattering length can be tuned by means of a Feshbach 
resonance and which has a /3@ value of 164a;, ^ r , where 
o-bohr denotes the Bohr radius. For this system, a scat- 
tering length of a s — 16.9/36 corresponds to 2770af, o / lr , 
a comparatively large value that can be realized experi- 
mentally in 85 Rb gases. As a point of reference, a density 
of 10~ 5 /3g~ 3 corresponds to a density of 1.53 x 10 13 cm -3 
for 85 Rb. 

The solid curve with positive energy in Fig. [2] shows 
the energy per particle Eg t o/N at unitarity, Eg,o/N w 

13.3h 2 n 2 J 3 jm [l2|. As seen in Fig. this unitary limit 
is approached by the energy per particle for the Bose 
gas with a s = 169/36 (filled circles connected by a dotted 
line). To illustrate this point, Fig. [3] shows the scaled 
average interparticle spacing r //3 6 as a function of the 
scaled density n_e/3f for a s = 169/36- This plot indicates 
that the unitary requirement, a s 3> r a 3> R, is met for 
values of UbPq larger than about 10~ 5 . Similarly, we find 



that the family of liquid curves converges to Eb,q/N w 

— 2A6h 2 n 2 J 3 /m (see Sec.lIVIfor details), plotted as a solid 
line in Fig. [21 when the inequalities a s 3> r a 06 are 
fullfilled. We note that the unitarity curve with negative 
energy is also approached, from above, for systems with 
large negative scattering lengths (not shown in Fig. \2§ ■ 
Aside from the proportionality constant, the power law 
relation for the liquid and gas branches at unitarity is 
the same. 

In addition to a Bose system interacting through the 
van der Waals potential, we consider a Bose system in- 
teracting through the square well potential with range 
R. For a given scattering length a s and density ng, the 
energy per particle Eb,q/N for these two two-body po- 
tentials is essentially identical for the densities shown in 
Fig. [2] This agreement emphasizes that the details of 
the two-body potential become negligible at low density, 
and in particular, that the behavior of the Bose gas in 
the unitary limit is governed by a single length scale, the 
average interparticle spacing r a . 

As discussed in SecHH the formal parallels between the 
LOCV method applied to bosons and fermions allows the 
energy per particle Ep,o/N for a two-component Fermi 
gas, Eq. (|17p . to be obtained straightforwardly from the 
energy per particle Eb,o/N of the Bose system. Fig- 
ured] shows the dimensionless energy (Ef,o/N)/Efg as 
a function of the dimensionless quantity l/(kpa s ) for the 
square well potential for np = 10~ 6 i? -3 . We find essen- 
tially identical results for the van der Waals potential. 
The crossover curve shown in Fig. 3] describes any dilute 
Fermi gas for which the range R of the two-body poten- 
tial is very small compared to the average interparticle 
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spacing r D . In converting the energies for the Bose sys- 
tem to those for the Fermi system, the gas branches of the 
Bose system (dotted lines in Figs. [2] and [3]) "turn into" 
the excited state of the Fermi gas (dotted line in Fig. Q|; 
the liquid branches of the Bose system with positive a s 
(dashed lines in Figs. [2] and [3]) "turn into" the part of 
the BEC-BCS crossover curve with positive a s (dashed 
line in Fig. [|]); and the liquid branches of the Bose system 
with negative a s (not shown in Figs. [2] and [3]) "turn into" 
the part of the BEC-BCS crossover curve with negative 
a s (dash-dotted line in Fig. 2]). 

To emphasize the connection between the Bose and 
Fermi systems further, let us consider the BEC side of 
the crossover curve. If l/(kpa s ) > 1, the fermion energy 
per particle Ep^/N is approximately given by Edimer/2, 
which indicates that the Fermi gas forms a molecular 
Bose gas. Similarly, the liquid branch of the Bose sys- 
tem with positive scattering length is made up of bosonic 
molecules as the density goes to zero. The formal analogy 
between the Bose and Fermi LOCV solutions also allows 
the energy per particle Ep,a/N at unitarity, i.e., in the 
l/(fcp \a s \) — > limit, to be calculated from the energies 
for large a s of the gas and liquid branches of the Bose 
system (solid lines in Fig. [2]). For the excited state of the 
Fermi gas we find Ep^/N ss 2>.92Epg, and for the low- 
est gas state we find Ep^/N m QAQEpo- These results 
agree with the LOCV calculations of Ref. [20(, which use 
an attractive cosh-potential and a (5-function potential. 
The value of OAQEpc is in good agreement with the en- 
ergy of 0A2Epc obtained by fixed-node diffusion Monte 
Carlo calculations I2U 01 • 



IV. BOSE AND FERMI SYSTEMS BEYOND 
s-WAVE AT UNITARITY 



This section investigates the unitary regime of Bose 
and Fermi systems interacting through higher angular 
momentum resonances. These higher angular momen- 
tum resonances are necessarily narrow @ , and we hence 
expect the energy-dependence of the generalized scatter- 
ing length ai(k) to be particularly important in under- 
standing the many-body physics of dilute atomic systems 
beyond s-wave. In the following we focus on the strongly- 
interacting limit. Figure [5] shows \/ai{k) as a function of 
the relative scattering energy E re i for the square- well po- 
tential with infinite zero-energy scattering length a; for 
three different angular momenta, I — (solid line), I = 1 
(dashed line), and I = 2 (dotted line). Figure [5] shows 
that the energy-dependence of ai(k) increases with in- 
creasing /. 

Our goal is to determine the energy per particle 
Eb,i/N for Bose systems with finite angular momentum 
I in the strongly-interacting regime. For s-wave interac- 
tions, the only relevant length scale at unitarity is the av- 
erage interparticle spacing r Q (see Sec. IIII[) . In this case, 
the energy per particle at unitarity can be estimated an- 
alytically by evaluating the LOCV equations subject to 



0.4- 



0.2 
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I 



-0.2- 
-0.4 



1=0 
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1=2 



-0.01 



-0.005 

E re ,/(h 2 /mR 2 ) 



0.005 



0.01 



FIG. 5: R/ai(E re i) as a function of the scaled relative scatter- 
ing energy E re i/(h 2 /mR 2 ) for the square well potential Vs-w 
with infinite zero-energy scattering length ai, i.e., l/ai = 0, 
for three different partial waves [1 = (solid line), I — 1 
(dashed line), and I = 2 (dotted line)]. 



the boundary condition implied by the zero-range s-wave 
pseudo-potential [l2j]. Unfortunatey, a similarly simple 
analysis that uses the boundary condition implied by the 
two-body zero-range pseudo-potential for higher partial 
waves fails. This combined with the following arguments 
suggests that Eba/N depends additionally on the range 
of the underlying two-body potential for finite I: i) The 
probability distribution of the two-body Z-wave bound 
state, I > 0, remains finite as a/ approaches infinity and 
depends on the interaction potential [1^, H(| ■ h) A de- 
scription of 1-wa.ve resonances (I > 0) that uses a coupled 
channel square well model depends on the range of the 
square well potential [4l[ . iii) The calculation of struc- 
tural expectation values of two-body systems with finite 
I within a zero-range pseudo-potential treatment requires 
a new length scale to be introduced (2f| . 

Motivated by these two-body arguments (see also 
Refs. [12, [H, E| for a treatment of p-wave interacting 
Fermi gases) we propose the following functional form for 
the energy per particle Eba/N of a 1-wa.ve Bose system 
at unitarity interacting through the square-well potential 
v sw with range R, 



E 



B.I 



N 



Ci 



mR x i/ 2 



2/3-xj/6 



(18) 



Here, C; denotes a dimensionless ^-dependent proportion- 
ality constant. The dimensionless parameter xi deter- 
mines the powers of the range R and the density ns, 
and ensures the correct units of the right hand side of 
Eq. ([HI). To test the validity of Eq. ([TgJ). we solve the 
LOCV equations, Eqs. {TTJ through flU}, for I = to 2 
for the one-component Bose system. Note that the one- 
component p-wave system is unphysical since it does not 
obey Bose symmetry; we nevertheless consider it here 
since its LOCV energy determines the energy of two- 
component p-wave Bose and Fermi systems (see below) . 
Figure [S] shows the energy per particle Eba/N for 
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FIG. 6: Scaled energy per particle (E B ,i/N)/(h 2 /mR 2 ) for a 
one-component Bose system for the energetically lowest-lying 
gas branch as a function of the scaled density ubR 3 obtained 
by solving the LOCV equations [Eqs. (fTTjl through (fl4|) ] for 
v aw for three different angular momenta I [I = (crosses), 
I = 1 (asterisks) and I = 2 (pluses)]. The depth Vo of v aw is 
adjusted so that l/ai(k) = 0. Solid, dotted and dashed lines 
show fits of the LOCV energies at low densities to Eq. (|18p 
for I = 0, I = 1 and I = 2 (see text for details). Note that 
the system with I = 1 is of theoretical interest but does not 
describe a physical system. 



lying gas branch. The density ranges used in the fit are 
chosen so that Eq. (fig)) describes the low-density or uni- 
versal regime accurately. Solid, dotted and dashed lines 
in Fig. [5] show the results of these fits for I = 0, 1 and 
2, respectively; in the low density regime, the lines agree 
well with the symbols thereby validating the functional 
form proposed in Eq. (|18p. 

We repeat the LOCV calculations for the energeti- 
cally highest-lying liquid branch of the one-component 
Bose system. By fitting the LOCV energies of the liquid 
branches for small densities to Eq. (118|) , we determine Xi 
and Ci [45| . We find the same xi as for the gas branch but 
different C\ than for the gas branch (the proportionality 
constants obtained for the liquid branch are denoted by 
C; L ; see Table [TJ> . Our values for xp, Cn and Cq agree 
with those reported in the literature [l2t l20l . l46j . 

Equation (|18[) can be rewritten in terms of the com- 
bined length £/, 



E 



B.l 



N 



(1/6-2/3) h 2 



where 



6 = r£- l '^R 1 / 4 



(19) 



(20) 



/ 


Xl 


c? 


c? 





0.00 


-2.46 


13.3 


1 


1.00 


-3.24 


9.22 


2 


2.00 


-3.30 


6.98 



TABLE I: Dimensionless parameters xi, and Cf for I = 
to 2 for a one-component Bose system obtained by fitting the 
LOCV energies Eb,i/N for small densities to the functional 
form given in Eq. (Tl8|) (see text for details). 



a one-component Bose system, obtained by solving the 
LOCV equations for the energetically lowest-lying gas 
branch, as a function of the density Ub for I = (crosses), 
I = 1 (asterisks), and I = 2 (pluses) for the square well 
potential, whose depth Vo is adjusted for each I so that 
the energy- dependent generalized scattering length ai(k) 
diverges, i.e., l/a;(fc) = 0. Setting aj(fc) to infinity en- 
sures that the Z-wave interacting Bose system is infinitely 
strongly interacting over the entire density regime shown 
in Fig. [5] Had we instead set the zero-energy scattering 
length ai to infinity, the system would, due to the strong 
energy-dependence of ai (k) [see Fig. [5], "effectively" in- 
teract through a finite scattering length. 

Table U summarizes the values for xi and Cj F, which we 
obtain by performing a fit of the LOCV energies Eb.i/N 
for the one-component Bose system for small densities to 
the functional form given in Eq. (| 18[) . In particular, we 

find xi — i, which implies that Eb,i/N varies as n 2 ^ 3 , 

n 1 ^ 2 and n 1 ^ 3 for I = 0, 1 and 2, respectively. Table U 
uses the superscript "G" to indicate that the proportion- 
ality constant is obtained for the energetically lowest- 



For the s-wave case, £j reduces to r D and the conver- 
gence to the unitary regime can be seen by plotting 
{Eb.o/N)/(K 2 /mr 2 ) as a function of a s /r D [12]. To 
investigate the convergence to the unitary regime for 
higher partial waves, Fig. [7] shows the energy per parti- 
cle Eb,i/N for the energetically lowest-lying gas branch 
as a function of ai(E re i)/£i for fixed energy-dependent 
scattering lengths a/(fc), i.e., for a s (k) — 10 10 i?, a p (k) 
I i) '" It. and a d (k) = 10 6 R [the different values of a t (k) 
are chosen for numerical reasons] . Figure [7] shows that 
the inequality 



\ai(E re i)\ >£ ; >i? 



(21) 



is fulfilled when (E b ,i /N) / (h 2 / mt;f) is constant. Note 
that this inequality is written in terms of the energy- 
dependent scattering length (see above). We find sim- 
ilar results for the liquid branches for I = to 2. For 
higher partial waves, we hence use the inequality given 
by Eq. (pH]) to define the unitary regime. In the uni- 
tary regime, the energy per particle Eb,i/N of the Bose 
system depends only on the combined length scale 
For s-wave interacting systems, we have £ s = r a and 
a s (k) m a s , and Eq. (|21[) reduces to the well known s- 
wave unitary condition, i.e., to |a s | 3> r Q 3> R. 

We now discuss those regions of Fig. [3 where the 
energy per particle (E b ,i / 'N) / '(ti 2 / 'm^ 2 ) for the one- 
component Bose system deviates from a constant. For 
sufficiently large densities, the characteristic length £/ be- 
comes of the order of the range R of the square well 
potential. In this "high density" regime the system 
exhibits non-universal behaviors. In Fig. e.g., the 
energy-dependent scattering length a^(k) equals 10 6 i?; 
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FIG. 7: Scaled energy per particle (E B ,i/ N) / (h 2 /mtf) for 
the energetically lowest-lying gas branch of /-wave interacting 
one-component Bose systems obtained by solving the LOCV 
equations [Eqs. pip through (I14|l ] for v sw as a function of 
ai (E re i)/£i ■ The depth Vq of v S m is adjusted so that a; (E re l) — 
10 10 R for I = and I = 1, and a t (E rel ) = 10 6 R for I = 2. 
Note that the system with I — 1 is of theoretical interest but 
does not describe a physical system. In the regime where the 
inequality R <C <C ai(E re {) is fulfilled, the scaled energy 
per particle is constant; this defines the unitary regime. 



correspondingly, ^ equals R when a<i(k)/£d = 10 6 - As 
a>d{k)/£,d approaches 10 6 from below, the system becomes 
non-universal, as indicated in Fig. [7] by the non-constant 
dependence of the scaled energy per particle on a^k)/^. 
On the left side of Fig. [TJ where ai(E re i)/£i becomes of 
order 1, the "low density" end of the unitary regime is 
reached. When a p (E re i)/^ p equals 10, e.g., the inter- 
particle spacing r a equals 10 2 a p (E re i), i.e., the system 
exhibits universal behavior even when the interparticle 
spacing is 100 times larger than the scattering length 
a p (E re i). This is in contrast to the s-wave case, where 
the universal regime requires \a s \ 3> r a . The different 
behavior of the higher partial wave systems compared to 
the s-wave system can be understood by realizing that 
is a combined length, which contains both the range R 
of the two-body potential and the average interparticle 
spacing r Q . For a given ai(E re i) / R, the first inequality 
in Eq. (|2"Tj) is thus satisfied for larger average interparti- 
cle spacings r a /R, or smaller scaled densities tibR 3 , as I 
increases from to 2. 

In addition to investigating /-wave Bose gases interact- 
ing through the square well potential v sw , we consider 
the van der Waals potential v v d w - For the energetically 
lowest-lying gas branch of the one-component "p-wave 
Bose" system we find the same results as for the square 
well potential if we replace R in Eqs. (fTS|) to (|2"U)) by 
(3q. We believe that the same replacement needs to be 
done for the liquid branch with I = 1 and for the liquid 
and gas branches of e?-wave interacting bosons, and that 
the scaling at unitarity derived above for the square well 
potential holds for a wide class of two-body potentials. 

Within the LOCV framework, the results obtained for 



/ 




ct 


c? 





0.00 


-1.55 


8.40 


1 


1.00 


-2.29 


6.52 


2 


2.00 


-2.62 


5.54 



TABLE II: Dimensionless parameters xi, C ; L and Cf for / = 
to 2 for a two-component Bose system (see text for details) . 



xi, Cf and for the one-component Bose systems can 
be applied readily to the corresponding two-component 
system by scaling the Bose density appropriately (see 
Sec. [II]). The resulting parameters xi, and Cf for 
the two-component Bose systems are summarized in Ta- 
ble nn 

The energy per particle Ep^/N for /-wave interact- 
ing two-component Fermi systems can be obtained from 
Eq. (fT7|) using the LOCV solutions for the liquid and 
gas branches discussed above for /-wave interacting one- 
component Bose systems. In the unitary limit, we find 



E 



F.I 



N 



A — n 



2/3 



B, 



h 2 n 2 J 3 - l/6 



R 1 / 2 



(22) 



where A = (3/10)(3tt 2 ) 2 / 3 w 2.87 and Bi = 
c .G,Zy 2 2/ 3 -i/6 ( the c G,l are giyen in Tablc rj| The first 

term on the right hand side of Eq. (|22|) equals Epc, and 
the second term, which is obtained from the LOCV solu- 
tions, equals A//2. The energy per particle Epj/N at uni- 
tarity is positive for all densities for Bi = Cp/2 2//3 ~'/ 6 . 
For Bi = C/ J /2 2 / 3_/ / 6 , however, the energy per particle 
Ef.i/N at unitarity is negative for / > for small den- 
sities, and goes through a minimum for larger densities. 
This implies that this branch is always mechanically un- 
stable in the dilute limit for / > 0. 

The LOCV treatment for fermions relies heavily on the 
product representation of the many-body wave function, 
Eq. (|15p . which in turn gives rise to the two terms on 
the right hand side of Eq. (|22|) . It is the competition of 
these two energy terms that leads to the energy minimum 
discussed in the previous paragraph. Future work needs 
to investigate whether the dependence of Epj/N on two 
length scales as implied by Eq. ([22]) is correct. In contrast 
to the LOCV method, mean-field treatments predict that 
the energy at unitarity is proportional to Ep G\kFr e ,i\, 
where r e j denotes a range parameter that characterizes 
the underlying two-body potential [42|, 0, |44| . 



V. CONCLUSION 

This paper investigates Bose and Fermi systems us- 
ing the LOCV method, which assumes that three- and 
higher-order correlations can be neglected and that the 
behaviors of the many-body system are governed by two- 
body correlations. This assumption allows the many- 
body problem to be reduced to an effective two-body 
problem. Besides the reduced numerical effort, this for- 
malism allows certain aspects of the many-body physics 
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to be interpreted from a two-body point of view. Further- 
more, it allows parallels between Bose and Fermi systems 
to be drawn. 

In agreement with previous studies, we find that the 
energy per particle "corrected" by the dimer binding en- 
ergy, i.e., Ef,o/N — Edimer/2, of dilute two-component 
s-wave Fermi gases in the whole crossover regime de- 
pends only on the s-wave scattering length and not on 
the details of the underlying two-body potential. Fur- 
thermore, at unitarity the energy per particle is given by 
E F ,o/N = 0A6E FG . This LOCV result is in good agree- 
ment with the energy per particle obtained from fixed- 
node diffusion Monte Carlo calculations, which predict 
E F .o/N = 0A2E FG pi, [H I3. This agreement may 
be partially due to the cancellation of higher-order cor- 
relations, and thus somewhat fortuitous. In contrast to 
Ref. [28| , we find that the liquid branch of bosonic he- 
lium does not terminate at low densities but exists down 
to zero density. 

For higher angular momentum interactions, we deter- 
mine the energy per particle of one- and two-component 
Bose systems with infinitely large scattering lengths. For 
these systems, we expect the LOCV formalism to pre- 
dict the dimensionless exponent xi, which determines 



the functional dependenc of Eb.i/N on the range R of 
the two-body potential and on the average interparticle 
spacing r Dl correctly. The values of the proportionality 
constants CP and C Z L , in contrast, may be less accurate. 
We use the LOCV energies to generalize the known uni- 
tary condition for s-wave interacting systems to systems 
with finite angular momentum. Since higher angular mo- 
mentum resonances are necessarily narrow, leading to a 
strong energy-dependence of the scattering strength, we 
define the universal regime using the energy-dependent 
scattering length a/(fc). In the unitary regime, the en- 
ergy per particle can be written in terms of the length 
£i, which is given by a geometric combination of r Q and 
R. The LOCV framework also allows a prediction for 
the energy per particle of two-component Fermi gases 
beyond s-wave to be made [see Eq. (|22l) ]. Although the 
functional form of the many-body wave function for two- 
component Fermi systems used in this work may not be 
the best choice, we speculate that the energy scales de- 
rived for strongly interacting Bose systems are also rele- 
vant to Fermi systems. 

This work was supported by the NSF through grant 
PHY-0555316. RMK gratefully acknowledges hospitality 
of the BEC Center, Trento. 
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